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Impact of symmetry breaking in networks of globally coupled oscillators 
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We analyze the consequences of symmetry breaking in the coupling in a network of globally 
coupled identical Stuart-Landau oscillators. We observe that symmetry breaking leads to increased 
disorderliness in the dynamical behavior of oscillatory states and consequently results in a rich variety 
of dynamical states. Depending on the strength of the nonisochronicity parameter, we find various 
dynamical states such as amplitude chimera, amplitude cluster, frequency chimera and frequency 
cluster states. In addition we also find disparate transition routes to recently observed chimera 
death state in the presence of symmetry breaking even with global coupling. We also analytically 
verify the chimera death region which corroborates the numerical results. The above results are 
compared with that of the symmetry preserving case as well. 
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I. INTRODUCTION 

The study of complex networks has attracted much 
attention in the fields of physics, chemistry, social sci¬ 
ences, etc. under various coupling topologies 0 - 0 - For 
instance collective phenomena in complex networks have 
been studied profoundly in the past decade which depend 
on both the peculiarity of the oscillators and the nature of 
the coupling and its strength. In particular much of the 
recent studies in nonlinear dynamics have been devoted 
to understand the dynamics of globally coupled oscilla¬ 
tors due its applicability in real world networks [ffl- 

Out of the various dynamical behaviors in identical 
nonlinear coupled systems, the existence of hybrid states 
(chimera states) of combining both synchronized and 
desynchronized behaviors jla4l8l | has received wide at¬ 
tention because of its connection with many real world 
applications such as uni-hemispheric sleep fl9j . neural 
networks (20] and so on. Moreover, chimera states have 
also been observed in maps j2lj , complex networks pH ], 
time discrete and continuous chaotic systems p3|. The 
emergence of chimera states in ensembles of oscillators 
has been studied theoretically 2dM321 for more than a 
decade and experimentally [331 - 4351 ] reported recently un¬ 
der nonlocal interactions. 

A few years ago Daido and Nakanishi in Ref. Q have 
observed a new interesting phenomenon called swing of 
synchronized states, while studying the inhomogeneity 
induced by the introduction of diffusive coupling in glob¬ 
ally coupled Stuart-Landau oscillators without symmetry 
breaking in the coupling. They found that the synchro¬ 
nized state which has been destabilized because of the in¬ 
crease in the coupling strength is found to be restabilized 
for further raise of it. The diffusion in globally coupled 
systems induces the synchronized state mediated by the 
so called cluster states. 

However recent studies suggest that coexistence be¬ 


havior of chimera states is also observed in globally cou¬ 
pled networks as in the case of systems with nonlocal 
coupling. Sethia and Sen have pointed out the emer¬ 
gence of amplitude mediated chimera states (AMC) even 
in globally coupled oscillators without symmetry break¬ 
ing in the coupling [H] in a system of Ginzburg-Landau 
oscillators as the specific example (note that the Stuart- 
Landau oscillators correspond to a special choice of a 
control parameter of this system). A question to ask is 
whether the amplitude mediated chimera state (AMC is 
also known as frequency chimera) can exist for global cou¬ 
pling with symmetry breaking. We analyze this question 
in the present work and report the existence of AMC in 
the case of symmetry breaking in the coupling also, which 
further leads to a rich variety of dynamical states than 
the case where the coupling is of symmetry preserving 
type. 

On the other hand, symmetry breaking instability in a 
network of coupled oscillators with nonlocal coupling also 
leads to the existence of stable inhomogeneous steady 
states. Interplay of nonlocality with symmetry breaking 
leads to new a dynamical state, namely chimera death, 
which was reported by Zakharova et al [37] recently. In 
the chimera death state, the oscillators in the network 
partition into two coexisting domains, where in one do¬ 
main neighboring nodes occupy the same branch of the 
inhomogeneous steady state (spatially coherent oscilla¬ 
tion death (OD)) whereas in the other domain neighbor¬ 
ing nodes are randomly distributed among the different 
branches of the inhomogeneous steady state (spatially 
incoherent OD). In the present study, we will also inves¬ 
tigate whether such chimera death states can exist in the 
globally coupled system with symmetry breaking as well. 
Indeed we show that such states do exist in this case. 

In this article, we are motivated to study the de¬ 
tailed dynamical behavior of the globally coupled Stuart- 
Landau oscillators in the presence of symmetry breaking 
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in the coupling and compare it with the case when the 
symmetry is preserved. We show that the influence of 
symmetry breaking in the coupling leads to increased dis¬ 
order in the nature of dynamical states and we observe 
that their regions also get widened. In [36j, Sethia and 
Sen observed that in the case of symmetry preserved cou¬ 
pling the synchronized state is mediated through the fre¬ 
quency chimera state. In addition to this we observe that 
the synchronized state is mediated through the ampli¬ 
tude chimera state (different from AMC) as well in both 
the cases of symmetry breaking and symmetry preserving 
couplings. We illustrate the above results with the help 
of characteristic measures such as standard deviation and 
strength of incoherence. Moreover, we present the re¬ 
sults for different transition routes to chimera death and 
also analytically verify the chimera death regions which 
nearly matches with the numerical results. 

This paper is organized as follows. In section II , we 
introduce the model of globally coupled Stuart-Landau 
oscillators with two different couplings that we have con¬ 
sidered for our simulation and present the results ob¬ 
tained from the analysis of the various dynamical states. 
In section III, we give the phase diagrams to illustrate 
the different transition routes. We summarize our find¬ 
ings in section IV. 


II. DYNAMICS OF GLOBALLY COUPLED 
STUART-LANDAU OSCILLATORS IN THE 
PRESENCE OF SYMMETRY PRESERVING AND 
SYMMETRY BREAKING COUPLINGS 

A. Model 

In order to exemplify our results, we consider an ar¬ 
ray of globally coupled identical Stuart-Landau oscilla¬ 
tors with two types of couplings: 

(i) Symmetry preserving coupling: 

Wj = Wj - (1 - ic)\wj\ 2 Wj + e(w - w-j ), ( 1 ) 

(ii) Symmetry broken coupling: 

Wj = Wj — (1 — ic)\wj\ 2 Wj + e(Re(w ) — Re(wj )), (2) 

where j = 1,2,3 N being the number of oscilla¬ 
tors, Wj = Xj + iyj, w = (1/N) SjLi w ji Re{w) = 

(1/N) x j is the mean field through the real part 
of the amplitude, e is the coupling constant and c is 
the nonisochronicity parameter. In our simulations, we 
choose the number of oscillators N to be equal to 100 
and in order to solve the equations © and © , we use 
the fourth order Runge-Kutta method with time step 
0.01. Note that Eq. © preserves the gauge symmetry 
Wj —► w'j = Wje‘ 9 , 9 € 1, or equivalently the rotational 
symmetry in_the (a :,y) plane, while it is broken in ©. 

In Ref. J|] Daido and Nakanishi have identified the 
phenomenon of swing in synchronized states in Eq. © 
which is mediated through cluster states. In our study we 
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FIG. 1: Standard deviation (a) of the system © with symme¬ 
try preserving coupling for different values of c (a) for c = 3.0, 
(b) c = 5.0, (c) c = 7.0, and (d-f) show their corresponding 
strength of incoherence S (red color line) and So (black color 
line). Regions-I, V correspond to synchronized states, regions- 
II, II' show the cluster states and the region-III represents 
chimera states. 


intensely study this phenomenon and clearly distinguish 
the dynamical regions in the following sections. We also 
compare the results with the system of oscillators with 
symmetry breaking in the coupling. It may also be noted 
that Eq. © is the special case of the Ginzburg-Landau 
equation 

Wj = Wj — (1 — ic)\wj\ 2 Wj + e(l + iC\)(w — Wj) (3) 

with the parameter chosen as Ci = 0. Occurrence of 
AMC in © in the (Ci,e) phase space is discussed in 
detail in Ref. [36} by Sethia and Sen. 


B. Quantification of various dynamical states 

We study the characteristic nature of the dynami¬ 
cal states including desynchronized, chimera/cluster and 
synchronized states in Eqs. © and © with the help of 
the quantitative measures such as standard deviation t 6] 
and strength of incoherence [38l [39 which are explained 
in the following subsections. 


1. CHARACTERIZATION WITH RESPECT TO 
STANDARD DEVIATION 

In this section, we use the measure standard deviation 
for the real part of the amplitudes as used by Daido and 
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FIG. 2: Standard deviation (cr) of the system © with sym¬ 
metry broken coupling for different values of c (a) for c = 2.5, 
(b) c = 4.0, (c) c = 7.0, and (d-f) show their correspond¬ 
ing strength of incoherence S (red color line) and So (black 
color line). Regions-I, ff correspond to synchronized states, 
regions-II, Iff show the cluster states and the region-III rep¬ 
resents chimera states. 


Nakanishi [§] defined by 

v={{\xj-x] | 2 ) 1/2 }, (4) 

where the bar denotes average over 1 < j < N and the 
angular bracket stands for the time average. It is clear 
that the standard deviation of the system that is cal¬ 
culated from d4j) is zero for the synchronized state and 
nonzero for the desynchronized state. 

Depending on the strength of the nonisochronicity pa¬ 
rameter c and coupling strength e, the system of oscilla¬ 
tors attains different dynamical states. In order to iden¬ 
tify the dynamical behavior of these states, to start with 
we make use of the above mentioned standard deviation 
measure for various strengths of c and e. With this inten¬ 
tion, we demonstrate the behavior of standard deviation 
a as a function of e in Figs, ©a-c) for symmetry preserved 
coupling (system ©) and in Figs, ©a-c) for symmetry 
broken coupling (system ©) for three different choices 
of c. 

(a) Symmetry preserved global coupling: 

In Fig. ©a), we fix the nonisochronicity parameter at 
c = 3.0 and one can observe that a takes nonzero val¬ 
ues for e < 0.275 which implies the desynchronization of 
oscillators in this region. For the values of e in the re¬ 
gion 0.275 < e < 0.675 the oscillators are synchronized 
where a decreases to zero value. As e increases a takes 
again nonzero values in the range 0.675 < e < 0.925 


which indicates desynchronization among the oscillators. 
By increasing e beyond 0.950, again a decreases to zero 
showing that the states correspond to a synchronized re¬ 
gion. Thus we can observe the recurrence of synchronized 
state for lower values of nonisochronicity parameter c. 
We can also observe the same behavior for higher values 
of nonisochronicity parameter. We illustrate this fact for 
c = 5.0 and c = 7.0 in Figs, [j] (b) and (c), respectively. 
This also confirms the existence of recurrence of synchro¬ 
nized states for higher values of the nonisochronicity pa¬ 
rameter as well. 

(b) Symmetry broken global coupling: 

To illustrate the dynamical behavior of the system © , we 
choose the nonisochronicity parameter as c = 2.5 which 
results in the dynamical states as depicted in Fig. ©a). 
For the values e < 0.27 the oscillators are in a synchro¬ 
nized state. By increasing the strength of coupling in¬ 
teraction (0.28 < e < 0.94), a decreases to zero so that 
the system of oscillators becomes synchronized. For the 
values of e in the region 0.95 < e < 2.48, er takes nonzero 
values implying that the oscillators in the system have 
again become desynchronized. By increasing e beyond 
2.48, we can observe synchronized states where er takes 
the value zero. Thus we can observe the recurrence of 
synchronized states for smaller values of nonisochronicity 
parameter c. We also confirm the existence of recurrence 
of synchronized states for higher values of nonisochronic¬ 
ity parameter c = 4.0 and c = 7.0 which is illustrated in 
Figs. [2] (b) and (c), respectively. Thus we can conclude 
that the system of oscillators that are coupled under sym¬ 
metry broken coupling follows the same transition as that 
of the system of oscillators with symmetry preserved cou¬ 
pling on the basis of the measure standard deviation er. 

However, we also note an important difference in the 
distribution of er as a function of e between the systems 
© and ©. While er varies smoothly as a function of the 
coupling strength e (Figs, ©a-c)) in the case of symme¬ 
try preserving global coupling corresponding to © , see 
also the work of Daido and Nakanishi [§], it varies quite 
randomly with symmetry broken case © as depicted in 
Figs. ©a-c). Thus we conclude that the presence of 
symmetry breaking in the system leads to an increase of 
disorder in the dynamical states. 


2. CHARACTERIZATION WITH RESPECT TO 
STRENGTH OF INCOHERENCE 

By using the above notion of standard deviation intro¬ 
duced in [6|, one cannot differentiate the chimera/cluster 
states and desynchronized states. In order to know the 
nature of dynamical states in more detail, we look at the 
strength of incoherence of the system a notion introduced 
recently by Gopal, Venkatesan and two of the present au¬ 
thors [38|] that will help us to detect interesting collective 
dynamical states such as the chimera state. For this pur¬ 
pose we introduce a transformation Zj = x 3 — Xj+\ 1381 ] . 
where j = 1,2, 3,..., IV. We divide the oscillators into M 
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FIG. 3: Snapshots for the variables Xj of the system Jl]) (with symmetry preserving coupling) and the corresponding frequencies 
/ for various dynamical states, (a)-(b) Amplitude cluster state for c = 5.0 and e = 0.60. (c)-(d) Amplitude chimera state for 
c = 5.0 and e = 0.75. (e)-(f) Frequency cluster state for c = 7.0 and e = 0.50. (g)-(h) Frequency chimera states for c = 7.0 and 

t = 1 . 0 . 


bins of equal length n = N/M and the local standard 
deviation <r(m) is defined as 


cr(m) = ((^ \zj - z j \ 2 ) 1/2 ) t ,m= 1,2, ...M. 

j=n(m— 1)+1 

(5) 

From this we can find the local standard deviation for 
every M bin of oscillators that helps to find the strength 
of incoherence [§8j through the expression 


S=l- 


l^m=l Sn 

M 


,s m = 0(S - a(m)), 


( 6 ) 


where 6 is the threshold value which is small and 0 is 
the Heaviside step function. When <t(to) is less than S , 
s m takes the value 1, otherwise it is 0. Thus the strength 
of incoherence measures the amount of spatial incoher¬ 
ence present in the system which is zero for the spatially 
coherent synchronized state. It has the maximum value, 


that is S = 1, for the completely incoherent desynchro¬ 
nized state and has intermediate values between 0 and 
1 for chimera states and cluster states. Further to dis¬ 
tinguish the chimera and cluster states, we make use of 
the quantitative measure So [39} which is the strength 
of incoherence before the removal of discontinuity points 
while S is calculated in this region after removal of such 
points. Here the value of Sq is the same as S for desyn¬ 
chronized, synchronized and chimera states, but for the 
cluster states Sq takes nonzero values between zero and 
one while S takes the value zero. For more details see 

Ref. [1} 

In order to understand the different dynamical states 
and their transitions more clearly, we analyze the 
strength of incoherence corresponding to the cases dis¬ 
cussed in terms of a in the above sub-section for the 
symmetry preserving case (Figs. QJd-f)) and symmetry 
breaking case (Figs. [2Kd-f)). 

(a) Symmetry preserving global coupling: 

Fig. [TJd) is plotted for the strength of incoherence (S 
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FIG. 4: Snapshots for the variables Xj of the system {2} (with 
symmetry broken coupling) and the corresponding frequen¬ 
cies / for various dynamical states, (a)-(b) Amplitude cluster 
state for c = 4.0 and e = 1.0. (c)-(d) Amplitude chimera 
state for c = 4.0 and e = 1.30. (e)-(f) Frequency cluster state 
for c = 7.0 and e = 0.6. (g)-(h) Frequency chimera state for 
c = 7.0 and e = 1.0. 


as well as So) with respect to the coupling strength e 
for c = 2.5. Initially the oscillators are desynchronized 
in the region e < 0.27 where S (and also So) takes the 
value unity. By increasing the coupling strength e, S (and 
also So) reaches the value zero where the oscillators are 
in synchronization (region-I). Further for the values of e 
between 0.675 and 0.98, So takes nonzero value (below 
one) while S takes the value zero (region-II) indicating 
the presence of cluster states. Interestingly, on further 
increasing e again we can observe the synchronized state 
where the value of (S, So) is zero (region-I'). Thus we can 
conclude that the synchronized state is mediated through 
cluster states for lower values of the nonisochronicity pa¬ 
rameter c. 

Further we address the question whether the synchro¬ 
nized state is mediated through cluster states for higher 
values of nonisochronicity parameters also. This is ana¬ 
lyzed by the behavior of S and So for the values c = 4.0 
and c = 7.0 which is depicted in Figs. [l](c) and[2f), re¬ 
spectively. When c = 4.0 (Fig. [2(e)), we can observe that 
for small values of e (e < 0.27) S and So take the value 
unity which shows the desynchronization of oscillators in 
this region and by increasing e to 0.28, the oscillators 
are synchronized where the value of S and So decreases 
to zero. We can observe the presence of cluster states 
(where S = 0 and 0 < So < 1, region-II) in the range 
0.96 < e < 1.23. Interestingly, further increase in e leads 
to the occurrence of chimera states where S (and also 
So) oscillates between zero and one (region-III). Again 


beyond e = 1.23 we can observe that the value of S = 0 
and So takes the value between zero and one confirm¬ 
ing the presence of cluster states (region-II'). Further 
increase of e leads to synchronization of oscillators with 
S = 0 (and also So = 0) in the region-I'. 

Next, a study of the strength of incoherence for c = 7.0 
confirms a similar transition behavior which is depicted 
in Fig. [2(f). Hence we conclude that the synchronized 
state is mediated through chimera states in addition to 
the cluster states for higher values of nonisochronicity 
parameter c. 

(b) Symmetry broken global coupling: 

Next, Figs. [2(d-f) depict the behavior of S (and So) for 
the system © with symmetry breaking in the coupling 
for three values of the nonisochronicity parameter. For 
c = 3.0, the oscillators follow the transition route (Fig. 
[2(d)) as desynchronization —> synchronization —> cluster 
states —»• synchronization. When c = 4.0 (Fig. [2(e)) and 
c = 7.0 (Fig. (2(f)) the transition route is represented as 
desynchronization —> synchronization —> cluster states —> 
chimera states —> cluster states —> synchronization. We 
also note here that the region of chimera states is much 
wider in the present case, Figs. (2d-f), compared to the 
symmetry preserving case as seen in Figs. 0d-f). Hence 
the characterization with respect to S also confirms that 
the symmetry breaking in the coupling leads to an in¬ 
crease in disorder in the dynamical states and also the 
regions are widened, while both the system 0 and © 
follow the same transition route. In addition symmetry 
breaking in the system 0 ultimately leads to a specific 
feature called chimera death which is demonstrated in 
the following sub-section. 


3. CHARACTERIZATION OF CLUSTER, CHIMERA, 
AND CHIMERA DEATH STATES 

To identify the nature of the cluster and chimera states 
(and also the chimera death state) in more detail, we 
look into the snapshots for the variables Xj and the cor¬ 
responding frequencies of the oscillators. The frequency 
of each of the oscillators is calculated from the expression 
fj = 2Trflj /AT, where j = 1,2, 3, ...N, flj is the number 
of maxima in the time series Xj of the j th oscillator with 
time interval AT. 

(a) Symmetry preserving coupling: 

At first we consider the system 0 for the three differ¬ 
ent values of the nonisochronicity parameter c (similar to 
the above discussed cases) which are illustrated in Figs. 
[2 Snapshots of the variables Xj and their corresponding 
frequencies fj (Fig. [2(a) and (b)) are plotted for c = 5.0 
and e = 0.60 confirming the presence of cluster states. 
In this case, the oscillators within the each distinct clus¬ 
ter have identical amplitudes and phases. However the 
frequencies of the oscillators in different groups are the 
same [Fig. [2 (b)]. This type of cluster states is designated 
as amplitude cluster states. By increasing e to 0.75, we 
can also observe the amplitude chimera state (the fluctu- 
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FIG. 5: (a)Space-time plot and (b) Snapshot of the variables Xj for N = 100 oscillators which confirm the existence of chimera 
death state for coupling strength e = 3.5 and nonisochronicity parameter c = 2.5. 
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FIG. 6: Phase diagram for the system JTJ (symmetry pre¬ 
served coupling). Green color shows the region of synchro¬ 
nized states, yellow color shows desynchronized state (DS) re¬ 
gion, blue color represents amplitude chimera state (AC), vi¬ 
olet color shows the frequency chimera state (FC) and brown 
color represents the amplitude cluster state (ACS) and fre¬ 
quency chimera state (FCS). Region ‘A’: DS —»ACS —¥ SY, 
region ‘B’: DS ->ACS-> AC ->ACS SY, and region ‘C’ : 
DS -s-FCS-> FC —>FCS SY. 


ations exist only in the amplitudes while the frequencies 
of all the oscillators are the same) which is shown in Figs. 
[3[c) and (d). 

Interestingly, a further increase in the nonisochronic¬ 
ity parameter leads to an increase in disorder in the fre¬ 
quency of the dynamical states as well. For c = 7.0, we 
can observe the frequency cluster states. This means that 
the oscillators within each cluster have identical ampli¬ 
tudes and phases and also the same frequencies. These 
are different for different groups as shown in Figs. de,f) 
for e = 0.50. This type of cluster states are designated 
as frequency cluster states. Figs. EKg) and (h) are the 
snapshot and frequencies of the oscillators corresponding 
to a frequency chimera state for e = 1.0. Here the fre¬ 


quencies of the oscillators in the coherent region are the 
same, while they are randomly distributed in the region 
corresponding to incoherent behavior [Fig. [21(h)]. Apart 
from the above described states, no other dynamical state 
(chimera death state) has been identified in this case of 
symmetry preserving coupling. 

(b) Symmetry broken coupling: 

Next we discuss the nature of cluster and chimera states 
for the system m in the case of symmetry broken cou¬ 
pling. The snapshots for the variables Xj and the fre¬ 
quencies of the oscillators which are shown in Figs. SKa) 
and (b) respectively show the presence of amplitude clus¬ 
ter states for c = 4.0 and e = 1.0. On further increasing e 
to e = 1.30, one finds the existence of amplitude chimera 
states (Figs. [He), (d)). For higher values of the non¬ 
isochronicity parameter (c = 7.0), we can observe the 
frequency cluster states for e = 0.6 (Figs. [4[e), (f)) and 
frequency chimera states for e = 1.0 (Figs. 0[g), (h)). 


Finally, while varying the coupling strength, the sym¬ 
metry breaking present in the coupling leads to chimera 
death state for larger values of e. It has the combined 
properties of chimera and oscillation death. The popula¬ 
tion of identical oscillators splits into two coexisting do¬ 
mains: (i) spatially coherent oscillation death (neighbor¬ 
ing oscillators populate in the same branch of inhomoge¬ 
neous steady state either as or x^) and (ii) spatially 
incoherent oscillation death (population of neighboring 
oscillators are completely random between and x^) 
as shown in Fig. [5] for c = 2.5 and e = 3.5, where x W 
and x^ represent the amplitudes of the oscillators cor¬ 
responding to the upper and lower branches of the in¬ 
homogeneous steady state, respectively. Here the total 
number of oscillators (100) splits into two equal groups 
of inhomogeneous steady state. Finally the diverse tran¬ 
sition routes to chimera death state is identified with a 
two parameter phase diagram in the next section. 
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FIG. 7: (a) Phase diagram for the globally coupled system © with symmetry breaking by varying the values of c and e for 
N = 100 oscillators. Green color shows the synchronized state, yellow color shows the state corresponding to desynchronized 
region (DS), blue color shows amplitude chimera state (AC), violet color shows the frequency chimera state (FC) and brown 
color shows the amplitude cluster states (ACS) and frequency chimera states (FCS) and red color shows the chimera death 
region. Region ‘A’: DS -aACS->- SY -A CD, region ‘B’: DS -aACS-> AC ^ACS -A SY -a CD, and region ‘C’ : DS -aFCS-> 
FC —>FCS —> SY —> CD. (b) shows the analytical plot for global coupling with symmetry breaking. Regions (I) and (IV) 
correspond to synchronized states, region (II) and (III) are for chimera and cluster states and the region (V) shows the 
chimera death states. The red color dot represents the Takens-Bogdanov point. 


III. DIFFERENT TRANSITIONS IN THE 
PARAMETRIC SPACE (e, c) 

We now present a comprehensive analysis of the differ¬ 
ent dynamical states and transition routes between them 
in the (e, c) two parameter space under both symmetry 
preserved and symmetry broken couplings. 

A. Global coupling under symmetry preservation 

In order to study the appearance of the swing of syn¬ 
chronized state in globally coupled oscillators intensely, 
we plotted the two parameter phase diagram in the para¬ 
metric space (e, c) in terms of the strength of incoherence 
S, as shown in Fig. [G] Different dynamical states for a 
given pair of values of e and c are identified by making 
use of the strength of incoherence (S). For this purpose 
we scan the ( e, c) plane with increments of 0.015 in e and 
0.045 in c. We choose random initial conditions between 
-1 to +1 for every choice of e and c values and observe 
the dynamics in the two parameter space in terms of S. 
For sufficiently small values of c, the system of oscillators 
attains a synchronized state directly from the desynchro¬ 
nized state. By increasing c, for c = 2.5 the synchronized 
state is mediated through amplitude cluster states and it 
follows the transition route as desynchronization —> syn¬ 
chronization —> amplitude cluster states —> synchroniza¬ 
tion. This route is illustrated in the region ‘A’. We can 
observe that the synchronized state is mediated through 
amplitude chimera states along with amplitude cluster 


states in the range 3.0 < c < 6.0 and the correspond¬ 
ing transition route is represented in the region ‘B’. This 
transition is the one reported by Daido and Nakanishi 
in They also identified that the cluster states exist 
as periodic and nonperiodic (quasi-periodic) desynchro¬ 
nized states. In addition we distinguish the cluster states 
as amplitude cluster states and amplitude chimera states 
and the transition route is represented as desynchroniza¬ 
tion —> synchronization —> amplitude cluster states —> 
amplitude chimera —> amplitude cluster states -A synchro¬ 
nization. Interestingly, on further increasing the value of 
c to 7.0, the synchronized states of the oscillators are me¬ 
diated through the frequency chimera states in addition 
to the frequency cluster states which is represented by 
the region ‘C’ as desynchronization —> synchronization 
-A frequency cluster states -A frequency chimera -A fre¬ 
quency cluster states -A synchronization. A summary of 
different transition routes is presented in Table I. 


B. Global coupling under symmetry breaking 

Next, for the case of globally coupled oscillators with 
symmetry broken coupling over a wide range of coupling 
strength, we plot the two parameter phase diagram in 
the parametric space (e, c) in terms of the strength of in¬ 
coherence S which is shown in Fig. [7] (a). One of the 
major distinguishing features in the present case is that 
the chimera death states occur for all values of the non- 
isochronicity parameter c for suitable values of e. Further 
the range of chimera states gets widened as a function of 








(7) 


e. Specifically we can observe that for sufficiently small 
values of c (c < 1.5) the system of oscillators attains 
chimera death state through a synchronized state from 
the desynchronized state. By increasing c to c = 2.5, 
the synchronized state gets mediated through the ampli¬ 
tude cluster state and attains the chimera death state. 
This route is represented in the region ‘A’ (Fig(7| (a)) 
and it follows the transition route as desynchronization 
—> synchronization —> amplitude cluster —> synchroniza¬ 
tion — > chimera death. On increasing the value of the 
nonisochronicity parameter to the range 2.5 < c < 5, we 
can observe the swing of synchronized state via ampli¬ 
tude chimera state (the oscillators in the coherent and 
incoherent regions having the same frequency) in addi¬ 
tion to the amplitude cluster state, but ultimately end¬ 
ing up in the chimera death state (beyond the region of e 
shown in Fig. 7 (a)). We can observe the appearance of 
amplitude chimera states in between the regions of the 
amplitude cluster states. Here it follows the transition 
route as desynchronization — > synchronization —> ampli¬ 
tude cluster -+ amplitude chimera —> amplitude cluster -+ 
synchronization —> chimera death (region ‘B’ in Fig. Q 
(a))). Interestingly, on further increasing c beyond c « 5, 
the nonisochronicity parameter induces a disorder in the 
frequencies of the dynamical states and causes the sys¬ 
tem to pass through a different set of dynamical states. It 
follows the transition route as desynchronization —> syn¬ 
chronization —> frequency cluster —> frequency chimera 
—> frequency cluster —>■ synchronization —> chimera death 
which is shown in the region ‘C’. 

In our study, we analyzed the dynamical system © 
over a wide range of nonisochronicity parameter and the 
main points are as follows. We identified the presence of 
amplitude chimera states. We have also observed that 
the synchronized states are mediated through the ampli¬ 
tude chimera states in addition to the cluster states in 
©• Further we identified that the chimera regions get 
widened due to the presence of symmetry breaking in the 
coupling. The chimera death states occur for sufficiently 
large e values for all values of c. We have also charac¬ 
terized the cluster states which appeared in the regions 
‘A’, ‘B’ and ‘C’ as amplitude and frequency clusters. The 
various transition routes are again summarized in Table 
I. 


C. Analytical results 

To understand the existence of the various states dis¬ 
cussed earlier analytically, we first consider the system 
© and then ©■ 

(i) Symmetry broken case: 

We start by assuming the chimera states correspond to 
the coexistence of synchronized and desynchronized iden¬ 
tical groups of oscillators so that the system © can be 
written as 

-jr =z s ~{ 1 - ic)\z s \ 2 z s + 


e{pRe(Z s ) + (1 - p)Re(Z d ) - Re(z s )), 

^ = z d - (1 - ic)\z d \ 2 z d + 
at 

e(pRe(Z s ) + (1 - p)Re(Z d ) - Re(z d )), (8) 


where s = 1,2,3.../, d = 1,2,3,... k , l + k = N, N being 
the number of oscillators, p = l/N , q = k/N , p + q = 1, 

z s = (VOELi- 2 ^ z d = (1/fc) Yll=\ Zd- Z s and z d are 
the states of the oscillators corresponding to synchro¬ 
nized and desynchronized states, respectively. Following 
refs. [3^, [40] , considering the p ~ 1 limit, we can see 
that the chimera state appears from the synchronized 
state by varying the coupling strength and we obtain 
z s = e lct and 


-jjjT = z d - (1 - ic)\z d \ 2 z d - e(Re(e lct ) - Re{z d )). 


The above equation has a solution z d = e lct that im¬ 
plies the completely synchronized manifold of the system 
which is always stable. In order to find the stability of the 
chimera/cluster state we apply a multi-time scale pertur¬ 
bation to the solution of z d so that z d = w(r)e lct , where 
w(t) represents the amplitude of the desynchronized os¬ 
cillators and t = to + T with to and r representing the fast 
and slow time scales, respectively. Thus the dynamics of 
each of the desynchronized oscillators is represented by 


dw 

dr 


(1 — ic)w 


(1 — ic)\w\ 2 w + ^(1 — w ) + 
\{1 -w))e~ 2 ' ct . 


(9) 


Averaging over the fast time scale to between 0 and 
we obtain 


-j- = (1 — ic)w — (1 — ic)\w\ 2 w + -(1 — w). (10) 
dr 2 


Eq. ([Hill has one stable fixed point at w = 1, which 
corresponds to complete synchronization manifold. The 
other two fixed points 


e — (1 + c ) ± V —(! + c 2 )e 2 + (1 + c 2 - e) 
1 ^ 2 , 31 =- 


2(1 + C 2 ) 


( 11 ) 


exist for 


1 + c 2 - J(1 + c 2 ) 3 

e < esN =- 9- 

c z 


( 12 ) 


From the linear stability analysis on equation m about 
the fixed points HID we get the Hopf bifurcation curve, 


2(-2 + V5 + 2c 2 + c 4 ) 

'» = -- ■ (13) 

The curves for saddle node and Hopf bifurcations are 
plotted by using Eqs. m and m which are shown 
by the red color solid and red color dotted lines in Fig. 
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mb). Saddle connection boundary is obtained by solv¬ 
ing m numerically and its boundary is given by the 
pink color line. Takens-Bogdanov point is denoted by 
etb = 2 (—1 + y/2) at |c| = 1. Thus, it turns out that for 
|c| > 1, the original system has a stable chimera/cluster 
state solution in the range esc < e < esN- In Fig. mb) 
regions I and IV correspond to synchronized states, re¬ 
gions II and III represent cluster and chimera states and 
the region V corresponds to the chimera death states as 
discussed below. 

On the other hand, chimera death represents the situ¬ 
ation where the total population is split into two groups 
of inhomogeneous steady states. Eqs © and © can be 
written as two populations of oscillators and is specified 
by 


= z h i - (1 - ic)\z h i\ 2 z h i + 

e(pRe(zhi) + (1 -p)Re(z h2 ) - Re{z hl )), (14) 

= z h2 - (1 - ic)\z h2 \ 2 z h2 + 

e(j>Re(z h i) + (1 - p)Re(z h2 ) - Re(z h2 )), (15) 

where Zhi and Zhi are the states of the oscillators corre¬ 
sponding to two groups of inhomogeneous steady states. 
Considering our numerical results, the above system 
has a trivial equilibrium point (aq, yi, x 2 , J/ 2 )=( 0 ,0,0,0) 
and a nontrivial (inhomogeneous) equilibrium point 
(x 1 ,y 1 ,x 2 ,y 2 )=(xl,yl,-xl,-yl) for p = 0.5 (as con¬ 
firmed numerically), where 


dz h i 

dt 

dz h2 

dt 


x 


* 

1 


2 V2cyl 
—6 + a 


(16) 


Vi = 


( , i „.\ / (-l+e)(e+a)+c 2 (-4+3e+a0 

(-e + a) -( 1 —c)2 e - 


2 y/2c 


(17) 


with a = y/4c 2 (—l + e) + e 2 . We find the above non¬ 
trivial fixed point is linearly stable for e > e c , where the 
critical value of e is given by 


e c = 


+ c ) 2 - 


(18) 


By using the above equation, we plotted the curve for 
the chimera death region (blue color line in Fig. ©a) and 
(b)) for two populations of oscillators. Fitting this curve 
with numerically plotted phase diagram, we can observe 
that analytical results closely match with the numerical 
results (see Fig. ©a)). 

(ii) Symmetry preserving case: 

Proceeding in a similar way for the symmetry pre¬ 
served global coupling system ©, Eqs. © and © can 
be rewritten as 


dz s 

dt 


= z s 


(1 - ic)\z s \ 2 z s + 


e(pz s + (1 - p)z d - z 8 ), (19) 

= z d - (1 - ic)\z d \ 2 z d + 
at 

e(pz s + (1 - p)z d - z d ). (20) 

As before, the slow scale dynamics of the each of the 
desynchronized oscillators can be represented as Q 

^ = (1 — ic)w — (1 — ic)|u;| 2 w; + e(l — w). (21) 
dr 

By replacing e by e/2 in the above equation, it can be 
reduced to Eq. m which corresponds to the case of 
symmetry breaking in the coupling. This implies that 
the dynamical regions in the system 0 get widened ap¬ 
proximately twice as that of the system ©. From m 
and (f20j) . looking for chimera death states, the two pop¬ 
ulation equations m and m can now be written as 


= Zhi - (1 - ic)\z h i\ 2 z h i + 

e{pz h i + (1 - p)z h2 - z h i), (22) 

= z h2 - (1 - ic)\z h2 \ 2 z h2 + 

e(pzhi + (1 - p)z h 2 - Zhi)- (23) 

We can check that the fixed points of the type m-m 
do not exist for Eqs. (E^ll and (E51) under symmetry pre¬ 
served coupling and so one cannot identify the existence 
of chimera death in the system ©■ 


dz h i 

dt 

dz h2 

dt 


IV. CONCLUSION 

In summary, we have investigated the common and dis¬ 
tinguishing features underlying the collective dynamics of 
globally coupled Stuart-Landau oscillators under two dif¬ 
ferent coupling schemes: (1). symmetry preserved cou¬ 
pling, (2). symmetry broken coupling. We have observed 
that the synchronized state is mediated through the am¬ 
plitude chimera states in addition to amplitude cluster 
states for lower values of the nonisochronicity parame¬ 
ter. For higher values of this parameter, the synchro¬ 
nized state gets mediated through the frequency chimera 
states in addition to frequency cluster states. Moreover, 
the presence of symmetry breaking in the coupling leads 
to increased disorder in the dynamical states and also 
leads to the widening of the various interesting dynam¬ 
ical regions. In addition we have also identified the ex¬ 
istence of chimera death states and diverse transitions 
routes to chimera death state in the case of symmetry 
broken coupling. 

Through a multi-time scale perturbation analysis, we 
have also analytically established the various regions of 
dynamical states including chimeras and chimera death 
states. We can thus conclude that symmetry breaking 
in global coupling leads to a rich variety of collective dy¬ 
namical states. 
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Appendix A: Transition routes between various 
collective dynamical states 


The different transition routes allowed by the coupled 
Stuart-Landau oscillators under the two different cou¬ 
pling schemes are summarized in table / for comparison 
purpose. 
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S. No 

Coupling schemes 

Transition routes 

1 

global coupling 
without symmetry breaking 

(i) DS SY 

(ii) DS —>ACS—>■ SY 

(iii) DS —>ACS—>■ AC -s-ACS -> SY 

(iv) DS ->FCS -5- FC—> FCS ->■ SY 

2 

global coupling 
with symmetry breaking 

(i) DS -> SY -s- CD 

(ii) DS ->ACS->- SY -* CD 

(iii) DS —>ACS—>■ AC ->ACS -> SY-)- CD 

(iv) DS —s-FCS -5- FC-> FCS —>• SY —>■ CD 


TABLE I: Different coupling schemes and corresponding transition routes. DS —» Desynchronized states, SY—» synchronized 
states, AC—> amplitude chimera states, FC —¥ frequency chimera states, CD —» chimera death states, ACS —> amplitude cluster 
states, FCS —» frequency cluster states. 




